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ABSTRACT 

We identify SDSSJ 121010.1+334722.9 as an eclipsing post-common-envelope binary, 
with an orbital period of P OID = 2.988 h, containing a very cool, low-mass, DAZ white 
dwarf and a low-mass main-sequence star of spectral type M5. A model atmosphere 



analysis of the metal absorption lines detected in the blue part of the optical spectrum, 
along with the GALEX near-ultraviolet flux, yields a white dwarf temperature of 
T effiW D = 6000 ± 200 K and a metallicity value of log[Z/H] = -2.0 ± 0.3. The 
Nai AA 8183.27,8194.81 absorption doublet is used to measure the radial velocity of the 
secondary star, K scc = 251.7 ± 2.0kms _1 and Fei absorption lines in the blue part of 
the spectrum provide the radial velocity of the white dwarf, A'wd = 95.3 ±2.1 km s _1 , 
yielding a mass ratio of q = 0.379 ± 0.009. Light curve model fitting, using the 
Markov Chain Monte Carlo (MCMC) method, gives the inclination angle as i = 
(79.05° - 79.36°) ± 0.15°, and the stellar masses as M WD = 0.415 ± 0.010 M Q and 
0^ ' M soc — 0.158 ± 0.006 Mq. Systematic uncertainties in the absolute calibration of the 

photometric data influence the determination of the stellar radii. The radius of the 
white dwarf is found to be i? W D = (0.0157 - 0.0161) ± 0.0003 R and the volume- 
averaged radius of the tidally distorted secondary is i? S oc, vol. aver. = (0.197 — 0.203) ± 
0.003 Rq. The white dwarf in SDSS1210 is a very strong He-core candidate. 
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1 INTRODUCTION 

Our understanding of stellar structure and evolution leads 
to the fundamental prediction that the masses and radii of 
stars obey certain mass-radius (M-R) relations. The cal- 
ibration and testing of the M-R relations requires accu- 
rate and model-independent measurements of stellar masses 
and radii, commonly achieved with eclipsing binaries (e.g. 
lAndersen|[l99ll : ISouthworth fc ClauserJl2007h . 

Among main-sequence (MS) stars, M-dwarfs of 
low mass (< 0.3M Q ), are the most ubiquitous. 



However, few ecl ipsing low-ma ss MS+MS binaries 
are known (e.g. | L6pez-Moralei 120071; iMorales et alj 



20091: ICakirh fc Ibanoglul I201C 



Irwin et al 



20101 ; 



Dimitrov fc Kiurkchieva 20101 and references therein) 



and have accurate measurements of their masses and radii, 
affecting the calibration of the low-mass end of the MS 
M-R relation. To further complicate matters, existing 
measurements consistently result in radii up to 15% larger 
and effective temperatur es 400 K or more below the value s 
predicted by theory (e.g. iRibasfeoOtj ; lLopez- Morales! 120071 '). 
This is not only the c ase for low-mass MS+MS binaries 
(|Bavless fc Oroszl 120061 ). but it is also present in field stars 
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ijBerger et al.ll2006l ; [Morales et alj|2008h and the host stars 
of transiting extra-solar planets ([Torres 20Q3). 

The situation is similar for white dwarfs, the most com- 
mon type of stellar remnant. Very few white dwarfs have 
mo del-independent meas urements of their masses and radii 
(see lParsons et al1l2010a| h and eclipsi ng WD+WD binaries 
have only recently been discovered ( Steinfadt et al.l l20ld : 
iParsons et ail l201ll ; iBrown et al] I2OIII) . Consequent ly, the 
finite temperature M-R relation of white dwarfs (e.g. IWoodl 
ll995l ; |Panei et al.l l2000h remai ns largely untested by obser- 
vations ( Provencal et al.l[l998l K 

An alternative approach leading to accurate mass and 
radius measurements for WDs and MS stars is the study of 
eclipsing WD+MS binaries. Until recently, the population of 
eclipsing WD+MS binaries had stagna ted with only seven 
systems known fsee lPvrzas et al1 [2009. for a list), a direct 
result of the small num ber of the entire WD+MS binaries 
sample (~ 30 systems; ISchreiber fc Gansic kel l2003h . 

However, in recent years, progress has been made 
thanks to the Sloan Digital Sky Survey (SDSS; 
lYork et all \200(j ). A dedicated search for WD+MS 
binaries co ntained in the spectros c opic SDSS Data 
Release 6 llAdelman-McCarthv et all 120081 ) and DR7 
l)Abazaiian et al] 20091') yielded mor e tha n 1600 sys- 
tems (e.g. iRebassa-Mansergas et all 120101 ). of which 
~ 1/3 are (sh ort-period) pos t -comm on-envelope bi- 
naries (PCEBs) ijSchreiber et all 120081 ). The majority 
of these PCEBs contain low-m ass, late-type M dwarfs 
l|Rebassa-Mansergas et ah! l2010h . while a large per- 
centage of the WD primaries are of low-mass as well 
|Rebassa-Mansergas et aill201lh . 

A significant fraction of eclipsing systems should exist 
among this sample of PCEBs. Identifying and studying these 
eclipsing systems will substantially increase the observa- 
tional constraints on the M-R relation of both WDs and MS 
stars. Therefore, we have begun the first dedicated search for 
eclipsing WD+MS binaries in t he SDSS, and 5 new systems 
have already been pub lished dNebot Gomez-Moran et alJ 
120091 ; iPvrzas et a.1.1 [20091 . but see also lDrake et al.ll2010l for 
a complementary sample). 

SDSSJ 121010.1+334722.9 (henceforth SDSS1210), the 
subject of this paper, is one of the new systems identified 
in this search. In what follows, we present our observations 
(Sec.[2|), determine the orbital period and ephemeris (Sec.[3|) 
and analyse the spectrum of the white dwarf (Sec. [4]). Radial 
velocity measurements (Sec. [5]) combined with light curve 
fitting (Sec. [6]) lead to the determination of the masses and 
radii of the binary components (Sec.[7|). We also explore the 
past and future evolution of the system (Sec. [8]). 



Table 1. SDSS coordinates and u, g, r, i, z magnitudes of the tar- 
get SDSS1210 and the comparison stars used in the analysis. We 
also provide the GALEX near-UV magnitude of SDSS1210. 



Star 


RA 


Dec 


u 


9 


r 


i 


2 


NUV 


T 


182.54221 


33.78969 


18.10 


16.94 


16.16 


14.92 


14.02 


20.821 


CI 


182.55470 


33.76832 


17.72 


15.98 


15.33 


15.11 


15.02 




C2 


182.54229 


33.73406 


19.95 


17.33 


16.01 


15.33 


14.94 




C3 


182.62616 


33.78141 


16.85 


15.80 


15.46 


15.34 


15.34 





Table 2. Log of the photometric and spectroscopic observations. 
For the LT observations, we also provide the number of one- hour 
observing blocks per night. 



Date 


Telescope 


Filter/Grating 


Exp. [s] 


Blocks Frames Eclipses 


2009 Apr 01 


LT 


V+R 


5 


1 


708 


1 


2009 Apr 02 


LT 


V+R 


5 


2 


1416 





2009 Apr 03 


LT 


V+R 


5 


2 


1416 


1 


2009 Apr 04 


LT 


V+R 


5 


2 


1416 


1 


2009 Apr 05 


LT 


V+R 


5 


3 


2124 


1 


2009 Apr 06 


LT 


V+R 


5 


1 


708 





2009 Apr 29 


WHT 


R600B/R1200R 


900 




1 




2009 May 02 


WHT 


R600B/R1200R 


900 




3 




2010 Apr 23 


WHT 


R600B/R1200R 


600 




1 




2010 May 18 


WHT 


R600B/R600R 


900 




12 




2011 Feb 06 


LT 


V+R 


5 


1 


720 


1 


2011 Mar 02 


LT 


V+R 


5 


1 


720 


1 


2011 Apr 02 


LT 


V+R 


5 


1 


720 


1 


2011 May 08 


LT 


V+R 


5 


1 


720 


1 


2011 Jul 03 


LT 


V+R 


5 


1 


720 


1 



with exposure times of 15-30 min taken over the course of 
three nights, we found large radial velocity variations that 
strongly suggested an orbital period of a few hours. We ob- 
tained time-series photometry of SDSS1210 with a 16-inch 
telescope equipped with an ST8-XME CCD camera, with 
the aim to measure the orbital period from the expected 
ellipsoidal modulation, and immediately detected a shallow 
eclipse in the light curve. Enticed by this discovery, we sched- 
uled SDSS1210 for additional high-time resolution photome- 
try, using RISE on the Liverpool Telescope (LT), with which 
a total of 9 eclipses were observed. 

Table[TJ lists the SDSS coordinates and magnitudes of 
SDSS1210 and the three comparison stars used in the 
analysis presented in this paper, while Tabled summarises 
our photometric and spectrosc opic observations. W e note 
that SDSS1210 has a GALEX (jMorrissev et alj|2007l ) near- 
ultraviolet (NUV) detection, but no far-ultraviolet (FUV) 
detection. 



TARGET INFROMATION, OBSERVATIONS 
AND REDUCTIONS 



SDSS 1210 was discovered by IRebassa-Mansergas et aD 
i|2010l ) as a WDMS binary dominated by the flux of a 
low-mass companion with a spectral type M5V, suggest- 
ing that the white dwarf must be very cool. Inspecting the 
Na 1 AA 8183.27,8194.81 doublet in the six SDSS sub-spectrcQ 

1 The sub-exposures that are co-added to produce one SDSS 
spectrum of a given object 



2.1 Photometry: LT/RISE 

Photometric observations were obtained with the robotic 
2.0 m Liverpool Telescope (LT) on La Palma, Canary Is- 
lands, using the high-sp eed frame-transfer CCD camera 
RISE dSteele et al.ll2004l) equip ped with a single wideband 
V+R filter (|Steele et all [20081 ). Observations were carried 
out in one hour blocks, using a 2 x 2 binning mode with ex- 
posure times of 5 seconds. 

The data were de-biased and flat-fielded in the 
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Table 3. Times of mid-cclipse (and their errors), O-C values (and 
their errors) and cycle number for the ephemeris of SDSS1210. 
Mid-eclipse times and errors are in MJD(BTDB), O-C values and 
errors are in seconds. 



0.02 
2454927.609 



Figure 1. Sample light curve of SDSS1210 obtained with a 5 sec. 
cadence using RISE on the LT on April 05, 2009. 



standard fashion within the LT reduction pipeline and 
aperture photometry wa s performed using SEXTRACTOR 
llBertin fc Arnoutsl Il996l ) in the manner described in 
iGansicke et all (|2004l ).~ 

A sample light curve is shown in Figure[T] The out-of- 
eclipse variation is ellipsoidal modulation, arising from the 
tidally deformed secondary. 



2.2 Spectroscopy: WHT/ISIS 

Time-resolved spectroscopy was carried out at the 4.2 m 
William Herschel Telescope (WHT) on La Palma, Canary 
Islands, equipped with the double-armed Intermediate dis- 
persion Spectrograph and Imaging System (ISIS). The spec- 
trograph was used with a 1" slit, and an 600 lines/mm grat- 
ing (R600B/R600R) on each of the blue and red arms, al- 
though a few spectra were obtained with a 1200 lines/mm 
grating on the red arm (R1200R). Both the EE VI 2 CCD 
on the blue arm and the REDPLUS CCD on the red arm 
were binned by three in the spatial direction and two in the 
spectral direction. This setup resulted in an average disper- 
sion of 0.88A per binned pixel over the wavelength range 
3643 — 5137A (blue arm) and 0.99A per binned pixel over 
the wavelength range 7691-9184A (red arm, R600R). From 
measurements of the full width at half maximum of arclines 
and strong skylines, we determine the resolution to be 1.4A. 

The spectra were reduced using the starlini^I p ackages 
KAPP A and FIGARO and then o ptimally extra cted l|Hornel 
Il986l ) using the PAMELif] code (|Marshl fl989l ). The wave- 
length scale was derived from Copper-Neon and Copper- 
Argon arc lamp exposures taken every hour during the ob- 
servations, which we interpolated to the middle of each of 
the science exposures. For the blue arm the calibration was 
determined from a 5th order polynomial fit to 25 lines, with 
a root mean square (RMS) of 0.029A. The red arm was also 
fitted with a 5th order polynomial, to 17 arclines. The RMS 
was 0.032A. 



2 Maintained and developed by the Joint Astronomy Centre and 
available from http://starlink.jach.hawaii.edu/starlink 

3 Available from htt p: //www. Warwick. ac.uk/go/trmarsh.] 



Mid-Eclipse [d] 


Error [d] 


O-C [s] 


Error [s] 


Cycle 


54923.0336744 


0.0000060 


-1 


1 





54925.0255324 


0.0000082 


1 


1 


16 


54926.1459281 


0.0000069 


-0 


1 


25 


54927.1418460 


0.0000087 


-0 


1 


33 


55599.1376175 


0.0000061 


3 


1 


5431 


55623.0396100 


0.0000056 


-1 


1 


5623 


55654.0375754 


0.0000081 


-0 


1 


5872 


55690.0151216 


0.0000063 


1 


1 


6161 


55745.9109933 


0.0000069 


-2 


1 


6610 



3 ORBITAL PERIOD AND EPHEMERIS 

We determined the orbital period and ephemeris of 
SDSS1210 through mid-eclipse timings. This was achieved 
as follows: 

Mid-eclipse times were measured by mirroring the ob- 
served eclipse profile around an estimate of the eclipse cen- 
tre and shifting the mirrored profile against the original un- 
til the best overlap was found. This method is particularly 
well-suited for the box-shaped eclipse profiles in (deeply) 
eclipsing PCEBs. 

An initial estimate of the cycle count was then obtained 
by fitting eclipse phases (0° bscrvcd _ ^j') -2 over a wide range 
of trial periods. Once an unambiguous cycle count was es- 
tablished, a linear fit, of the form T — To + P or b x E, was 
performed to the times of mid-eclipse versus cycle count, 
yielding a preliminary orbital ephemeris. 

Subsequently, we phase-folded our data set using this 
preliminary ephemeris and proceeded with the light curve 
model fitting (see Sec. [6}. Having an accurate model at hand, 
we re-fitted each light curve individually. This provides a 
robust estimate of the error on the mid-eclipse time, as our 
code includes the time of mid-eclipse To as a free parameter. 

Repeating the cycle count determination and the linear 
ephemeris fitting, as described above, we obtain the follow- 
ing ephemeris for SDSS1210, 

MJD(BTDB) = 54923.033 686(6) + 0.124 489 764(1) E (1) 

calculated on a Modified Julian Date-timescale and cor- 
rected to the solar system barycentre, with the num- 
bers in parentheses indicating the error on the last 
digit. Thus, SDSS1210 has an orbital period of P or b = 
2.987754 336(24) h. The mid-eclipse times, the observed mi- 
nus calculated values (O-C) and their respective errors are 
given in Tabled Given the short baseline, there is as yet 
no evidence for pe riod changes which a re frequently seen in 
such binaries (e.g. iParsons et al.ll2010bO . 



4 SPECTROSCOPIC ANALYSIS 

Whereas the SDSS spectrum of SDSS1210 remained in- 
conclusive with respect to the nature of the white dwarf 
IjRebassa-Mansergas et al.l 120101 ) . our blue-arm WHT spec- 
troscopy immediately revealed a host of narrow metal lines 
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4100 4200 4300 4400 

Wavelength [A] 

Figure 2. The normalised average WHT spectrum in the white 
dwarf restframe, along with line identifications for absorption 
lines originating in the white dwarf photosphere. 



that exhibit radial velocity variations anti-phased with re- 
spect to those of the M-dwarf. The WHT spectra obtained 
in May 2010, averaged in the white dwarf restframe and 
continuum-normalised, are shown in Fig.[2]and illustrate the 
wealth of absorption lines from Mg, Al, Si, Ca, Mn, and Fc. 
Similar metal lines have been detected in t he optical spectra 
of a few other coo l PCEBs, e.g. RRC ae (jZuckerman et all 
l200St ) or LTT560 jTappert et alj|2007ft . and indicate accre- 
tion of mass via a wind from the M-dwarf. 

We have analysed the blue WHT spectra using 
hydrogen-dominated but metal-polluted (DAZ) spectra cal- 
culated with the stellar atmosphere code described by 
iKoesterl (|2010T ). We fixed the surface gravity to logg = 7.70, 
as determined from the fits to the LT light curve (Sect. [5]). 
The model grid covered effective temperatures 5400 K ^ 
Teff ,wd ^ 7400 K in steps of 200 K and metal and He abun- 
dances of log[Z/H] = -3.0, -2.3, -2.0, -1.3, -1.0, with all 
relevant elements up to zinc included, and fixed their relative 
abundances ratios to the respective solar values. We then fit- 
ted the model spectra to the average WHT spectrum in the 
range 3645-3930 A, where the contribution of the M-dwarf is 
entirely negligible. A good fit is found for T e fi ,wd — 6000 K 
and metal abundances at ~ 0.01 their solar values, how- 
ever, the effective temperature and the metal abundances 
are strongly correlated (Fig.[3]). 

This degeneracy is lifted by including the GALEX de- 
tection of SDSS1210, as the predicted near-UV flux is a 
strong function of the effective temperature. The uncer- 
tainty in the absolute flux calibrations of our WHT spec- 
tra and the GALEX observations introduces a small sys- 
tematic uncertainty on the final result, and we settle for 
Tcff.wD = 6000 ± 200 K and log[Z/H] = -2.0 ± 0.3 . In- 
dependently, the weakness of the Balmer lines in the WHT 
spectrum also requires that T c ff ,wd < 6400 K. The spectral 
modelling of SDSS1210 is illustrated on Fig.g] 

Adopting the white dwarf radius from the light curve 
fit (Sect .[6] and [7]), 7?wd = 0.0159 R©, the flux-scaling fac- 
tor of the best-fit spectral model implies a distance of d ~ 
50 ± 5pc, wh ich is in good agreement with d ~ 66 ± 34 pc 
estimated by iRebassa-Mansergas et ail |2010l ) from fitting 
the M-dwarf. 




Figure 3. Results of model spectra fitting to the average WHT 
spectrum. The single, big point indicates the best fit solution. The 
contours indicate the regions where the x 2 of the fit is within 1, 2 
and 3 a (dotted, short-dashed and long-dashed lines respectively) 
of the minimum (single point). 



The detection of metals in the photosphere of the 
white dwarf allows an estimate of the accret ion rate (e.g. 
iDupuis et all Il993l ; iKoester fc Wilkenl |2006t ), as long as 
the system is in accretion-diffusion equilibrium. In cool, 
hydrogen-rich atmospheres, such as the one in SDSS1210, 
the diffusion time scales of the different metals detected in 
the WHT spectrum vary by a factor of ~ 2 for a given 
temperature, and are, for T b h,wd = 6000 K, in the range 
30000-60000 yearfl It is plausible to assume that the aver- 
age accretion rate over the diffusion time scales involved is 
constant, as the binary configuration (separation of the two 
stars, Roche-lobe filling factor of the companion) changes on 
much longer time scales. Summing up the mass fluxes at the 
bottom of the convective envelope, and taking into account 
the uncertainties in T ff and the metal abundances, gives 
M ~ (5 ± 2) x 10~ 15 M yr" 1 . There are now three PCEBs 
with similar stellar components that have measured accre- 
tion rates, RRCae (M ~ 4 x 10" 16 M s yr" 1 ; [DebesI 
LTT560 (M ~ 5 x 10" 15 M m yr" 1 ; iTappert et al l 
and SDSS1210 (M ~ 5 X 10" 15 M 



2006 



2011) 



©yr )■ 

Whereas SDSS1210 and LTT560 have similar orbital 
periods, the period of RR Cae is roughly twice as long, sug- 
gesting that the efficiency of wind-accretion decreases as the 
binary separation and Roche-lobe size of the companion in- 
crease, as is expected. A more systematic analysis of the 
wind-loss rates of M-dwarfs and the efficiency of wind ac- 
cretion in close binaries would be desirable, but will require 
a much larger sample of systems. 



4 For completeness, we note that because we have adopted solar 
abundance ratios for the metals, these small differences in dif- 
fusion time scales imply slightly non-solar ratios in the accreted 
material. In principle, the individual metal-to-metal ratios can 
be determined from the observed spectrum of the white dwarf, 
and hence allow to infer the abundances of the companion star, 
however, this requires data with substantially higher spectral res- 
olution to resolve the line blends. 
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Figure 4. Spectral modelling of SDSS1210. Main panel: the SDSS spectrum (black) and the GALEX NUV flux (black point), along 
with the best-fit white dwarf model (red*, T e g we> = 6000 K, logg = 7.70, log[Z/H] = —2.0) and the best-fit M-dwarf template for 
the companion (red*, spectral type M5). Inset (a): The sum of the white dwarf model and M-dwarf template provide a good match to 
the blue end of the SDSS spectrum (black), with the low flux of the M-dwarf template being the dominant limitation. Inset (b): Best-fit 
white dwarf model (red*) and the average WHT spectrum for A < 3930 A, where the M-dwarf contributes practically nothing to the 
observed flux. Inset (c): Best-fit white dwarf model (red*) and the average WHT spectrum (black) illustrating the weakness of the H/9 
and H7 lines of the white dwarf. Increasing the temperature very rapidly results in Balmer lines and/or a NUV flux that are inconsistent 
with the observations. *The coloured Figure is available in the online version only. 



5 THE SPECTROSCOPIC ORBIT 

Radial velocities of the binary com- 
ponents have been measured from the 
Fe 1 AA 4045.813,4063.594,4071.737,4132.058,4143.869 
absorption lines for the white dwarf and the 
Na 1 AA 8183.27,8194.81 absorption doublet for the sec- 
ondary star. 

The Fe 1 lines were simultaneously fitted with a second- 
order polynomial plus five Gaussians of common width and 
a separation fixed to the corresponding laboratory values. A 
sine fit to the radial velocities, phase-folded using the orbital 
ephemeris (Equation[T} yields Kwn = 95.3 ± 2.1 kms^ 1 
and 7wd = 24.2 ± 1.4 km s" 1 . 

The Nai doublet was fitted with a second-order poly- 
nomial plus two Gaussians of common width and a sepa- 
ration fixed to the corresponding laboratory value. A sine 
fit to the radial velocities, phase-folded using the orbital 



ephemeris yields K BCC — 251.7 ± 2.0 km s 1 and 7 SCC = 
12.2 ± 0.9 kms" 1 . 

Figure[5] shows the measured radial velocities phase- 
folded on the orbital period and the corresponding sine-fits. 

Knowledge of both radial velocities allows us to obtain 
the mass ratio q of the binary, namely q — Kwu/K scc — 
0.379 ± 0.009. We tentatively interpret the difference be- 
tween 7wd and 7 SCC as the gravitational redshift of the 
white dwarf zwD.spcc, which yields zwD.spcc = 11.9 ± 
1.7kms _1 (see also Sec.0. 



6 LIGHT CURVE MODELLING 

To obtain the stellar parameters of the binary components, 
light curve models were fitte d to the data using lcurve 
(see ICopperwheat et all [2010 for a description, as well as 
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Orbital Phase 

Figure 5. Phase-folded radial velocity curves of the secondary 
star (filled circles) and the white dwarf (open circles), with their 
respective errors. Also shown are the sine fits to the velocities of 
both components. A full cycle is repeated for clarity. 



iPvrzas et aill2009l ; iParsons et ai1l2010al ; iParsons et ai]|201ll 
for further applications). 



6.1 Code input 

The code computes a model based on input system parame- 
ters supplied by the user. The physical parameters defining 
the models are (i) the mass ratio q = M sec / A'lwu, (h) the 
binary inclination i, (iii) the stellar radii scaled by the bi- 
nary separation rwD = Rwu/a and r scc = R scc /a, (iv) 
the unirradiated stellar temperatures Teff.wD and T soc , (v) 
the sum of the unprojected stellar orbital speeds Vs = 
(-Kwd + K scc ) /sini, (vi) the time of mid-eclipse of the 
white dwarf To, (vii) limb- and gravity darkening coefficients 
and (viii) the distance d. The code accounts for the distance 
simply as a scaling factor that can be calculated very rapidly 
for any given model, and so it does not enter the optimisa- 
tion process. All other parameters can be allowed to vary 
during the fit. 



6.3 Minimisation 

Initial minimisation is achieved using the do wnhill- 
SIMP LEX and levenberg-marquardt methods (|Pressl 
|2002| ). while the Markov Chain Monte Carlo (MCMC) 
method (|Press et al.ll2007l ) was used t o determine the distri- 
butions of our model parameters (e.g. Ford 2006, and refer- 
ences therein). 

The MCMC method involves making random jumps 
in the model parameters, with new models being accepted 
or rejected according to their probability computed as a 
Bayesian posterior probability (the probability of the model 
parameters, 6, given the data, D, P (8\D)). P (8\D) is driven 
by a combination of \ 2 and a prior probability, P{9), that 
is based on previous knowledge of the model parameters. 

In our case, the prior probabilities for most parameters 
are assumed to be uniform. The photometric data provide 
constraints for the radii and inclination angle, however, the 
photometry alone cannot constrain the masses, as the light 
curve itself is only weakly depended on q. To alleviate this, 
we can use our knowledge of -Kwd and K sec . At each jump, 
the model values -Kw D and K!£ c are calculated through q, 
i and Vs. P (0) is then evaluated on the basis of the ob- 
served Kwo and K scc , assuming a Gaussian prior probabil- 
ity P(fi,a 2 ), with fi and a corresponding to the measured 
values and errors of JfwD and K sac . 

A crucial practical consideration of MCMC is the num- 
ber of steps required to fairly sample the parameter space, 
which is largely determined by how closely the distribution 
of parameter jumps matches the true distribution. We there- 
fore built up an estimate of the correct distribution starting 
from uncorrelated jumps in the parameters, after which we 
computed the covariance matrix from the resultant chain of 
parameter values. The covariance matrix was then used to 
define a multivariate normal distribution that was used to 
make the jumps for the next chain. At each stage the actual 
size of the jumps was scaled by a s ingle factor set to d eliver 
a model acceptance rate of sa 25% l|Roberts et al.|[l997l ) . Af- 
ter 3 such cycles, the covariance matrix showed only small 
changes, and at this point we carried out the long "pro- 
duction runs" during which the covariance and scale factor 
which define the parameter jumps were held fixed. 



6.2 Free and fixed parameters 

During the minimisation, we kept T^wd fixed at T e ff,WD = 
6000 K. The gravity darkening of the secondary was also kept 
fixed at 0.08 (the usual value for a convective atmosphere). 
Limb darkening coefficients were also held fixed. For the 
white dwarf we calculated quadratic limb darkening coef- 
ficients from a white dwarf model with T e ff ,wd = 6000 K 
and logg = 7.70, folded through the RISE filter profile. 
The corresponding values were found to be a = 0.174 and 
b = 0.421 for = l-a(l - fi) - 6(1 - ^) 2 , with 

H being the cosine of the angle between the line of sight and 
the sur face normal. For the seco ndary star we used the Ta- 
bles of IClaret fc Bloemenl l|201ll ). We interpolated between 
the values of V and R for a T = 3000 K and log g — 5 star, 
to obtain quadratic limb darkening coefficients a' = 0.62 
and b' = 0.273. All other parameters were allowed to vary. 



6.4 Stellar parameters 

Using the following set of equations, the stellar and binary 
parameters are obtained directly from the posterior distribu- 
tion of the model parameters, as outputed from the MCMC 
minimisation. 

The binary separation is obtained from the model parameter 
Vs through 

Porh 

a ~ "2V 



Vs 



(2) 



The white dwarf and secondary masses are obtained from 
the model parameters q and Vs as 



Mwd = 



and 



-Porb 

2ttG 1 



V S 3 



Msec = - — — — Vs 
2 ty G 1 + q 



(3) 



(4) 
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The stellar radii are directly obtained from the model pa- 
rameters rwD and r scc and Eq.[2]and the surface gravity of 
the white dwarf is of course given by 



logg = log 



/GMwd\ 



(5) 



6.5 Intrinsic data uncertainties 

The acquisition of high-precision absolute photometry on 
the LT in service mode is somewhat difficult to achieve. 
Each observing block individually covered only a third of 
the orbital phase and the blocks were obtained over many 
nights, under varying conditions (seeing, sky brightness, ex- 
tinction, airmass). The data are sensitive to changes in con- 
ditions, as they have been obtained through the very broad 
and non-standard V+R filter of RISE. In the absence of a 
flux standard, the photometry cannot be calibrated in ab- 
solute terms. When phase-folding the LT data, significant 
scatter is found at orbital phases where individual observing 
blocks with discrepant calibrations contribute. This affects 
both the shape of the eclipse, mainly the steepness of the 
WD ingress/egress and, to a lesser extend, the eclipse dura- 
tion, and the out-of-eclipse variation, i.e. the profile of the 
ellipsoidal modulation. As a result, there is an unavoidable 
systematic uncertainty in the photometric accuracy of our 
data, which will influence the determination of the stellar 
parameters. 

To gauge the effect of the systematic uncertainties we 
worked in the following fashion: each observing block has 
been reduced thrice, each time using one of the three com- 
parison stars reported in TableQ] CI has a g — r colour index 
comparable to SDSS1210, C2 is fairly red, while C3 is fairly 
blue. The data of each reduction were then phase-folded to- 
gether and two light curves were produced: one containing 
all the photometric points and one where (2-3) observing 
blocks with an obviously large intrinsic scattering were omit- 
ted. Thus, we ended up with six phase-folded light curves. 
A dedicated MCMC optimisation was calculated for each 
light curve. We will use the following notation when refer- 
ing to these chains: CIA denotes a light curve produced 
with comparison star CI and all data points, C2E denotes 
a light curve produced with comparison star C2 excluding 
observing blocks, and so on. 



7 RESULTS 

The results of the six MCMC processes are summarised in 
Tablef?] The quoted values and errors are purely of statisti- 
cal nature and represent the mean and RMS of the posterior 
distribution of each parameter. The radius of the secondary, 
as determined by r sec and a, is measured along the line 
connecting the centres of the two stars and, due to the tidal 
distortion, its value is larger than the average radius. There- 
fore, on TableUwe also report the more representative value 
of the volume-averaged radius. 

To illustrate the achieved quality of the fits, we plot 
models CIA and C1E in Figure[6] While the overall quality 
of the fit is very satisfactory, the model seems to slightly 
overpredict the flux at the "wings" of the ellipsoidal mod- 
ulation profile (phases ~ 0.05 — 0.15 and ~ 0.85 — 0.95). 




0-4 0-6 0-8 1 1.8 1-4 

White dwarf mass [M Q ] 

Figure 7. M ass-radius plot f o r wh it e dwarfs. Black points 
are data from IProvencal et al ] l|l998l l. IProvencal et al] fcoOih 
and ICasewell et al. I teoosh . The dotted line is the zero- 
temperature mass-radius re lation of Egglcton as quoted in 
IVerbunt fc Rappaport] Jl98Sft . The dashed line, marked as (He, 6) 
is a M-R relation for a T e g wd = 6000 K, He-core WD, with a hy- 
drogen lay er of M(H)/Mwp = 3 X 10~ 4 , interpo lated from the 
model s of lAlthaus fc Benvenutol il997l ). NNSer j Parsons et al] 
2010a) is marked, along with the track for T G ff wd = 60000 K, 
C/O-core WD, M(H)/M WD = 10" 4 (long dash-dot line), in- 
dicating the accuracy achieved in eclipsing PCEBs. The results 
of the six chains for SDSS1210 are plotted in red (online version 
only). Inset panel: zoom-in on the values of SDSS1210. The points 
are CIA: open circle; C1E: filled circle; C2A: open square; C2E: 
filled square; C3A: open triangle; C3E: filled triangle. 



This discrepancy could be data related, due to the intrin- 
sic scattering of points; system related, e.g. due to the pres- 
ence of starspots affecting the modulation; model related, as 
the treatment of stellar temperatures is based on blackbody 
spectra, for one specific wavelength; or due to a combination 
of these factors. 

With regard to the binary and stellar parameters, the 
MCMC results indicate the following: as expected for a de- 
tached system, the light curves depend very weakly on q and 
its value is well constrained by the radial velocitites. All six 
chains give inclination angle values just above 79°, consis- 
tent with each other within the errors. There is a slight shift 
upwards when excluding blocks from the phase-folded light 
curve. 

The tight spectroscopic constraints, mean that the com- 
ponent masses are largely independent of the model/data 
set used. Thus, the white dwarf in SDSS1210 has a mass of 
Mwd = 0.415 ± 0.010 Mq and the secondary star a mass 
of M scc = 0.158 ± 0.006 M . 

The quantity most seriously affected by systematics is 
the white dwarf radius. This is especially evident when con- 
sidering models C3A and C3E. However, such a discrepancy 
is expected, since C3 is considerably bluer than SDSS1210 
and is more susceptible to airmass/colour effects, leading to 
large intrinsic scattering. The values for Awd as obtained 
from CIA, C1E, C2A and C2E are consistent within their 
errors, indicating a systematic uncertainty comparable to 
the statistical one. This is illustrated in Figure[7] 

The secondary star radius is affected in a similar, albeit 
less pronounced, way. All six models lead to values broadly 
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Figure 6. Light curve fitting results for models CIA (left) and C1E (right). In each of the six panels we plot the phase-folded light curve 
with the model superimposed (top trace), the residuals of the fit (middle trace, offscted from for clarity) and a binned version of the 
residuals (bottom trace). Shown are the entire light curve (top panels), a zoom around the eclipse (middle panels) and the out-of-eclipse 
ellipsoidal modulation (bottom panels). 

Table 4. Stellar and binary parameters obtained from MCMC optimisation. The quoted values and errors are the mean and RMS of 
the posterior distribution of each parameter. The chains represent light curves created using comparison stars CI, C2 or C3 and either 
including all (A) observing blocks or excluding (E) those with obviously large scattering. See text for details. 



Parameter 



CIA 



C1E 



C2A 



C2E 



C3A 



C3E 



M WD [M ] 
%D [R©] 
WD logg 
A/ scc [M©] 

Usee [R©] 
-^scc, vol. aver. 

T scc [K] 



0.380 ± 0.010 

79.05 ± 0.15 

0.415 ± 0.010 

0.0157 ± 0.0003 

7.664 ± 0.015 

0.158 ± 0.006 

0.217 ± 0.003 

0.202 ± 0.003 

~ 2530 



0.380 ± 0.010 
79.28 ± 0.15 
0.414 ± 0.010 
0.0159 ± 0.0003 
7.652 ± 0.016 
0.157 ± 0.006 
0.212 ± 0.003 
0.199 ± 0.003 
~ 2550 



0.381 
79.03 
0.415 
0.0161 
7.641 
0.158 
0.217 
0.203 



Binary separation [R©] 0.871 ± 0.008 0.870 ± 0.008 0.871 



± 0.010 
± 0.15 
± 0.010 
± 0.0003 
± 0.015 
± 0.006 
± 0.003 
± 0.003 
2530 
± 0.008 



0.380 
79.13 
0.415 
0.0159 
7.649 
0.158 
0.215 
0.201 

0.871 



± 0.010 
± 0.15 



0.010 
0.0003 
0.017 
0.006 
0.003 
0.003 
2550 
± 0.008 



0.378 
79.36 
0.414 
0.0138 
7.773 
0.156 
0.210 
0.197 

0.869 



± 0.010 
± 0.18 
± 0.010 
± 0.0003 
± 0.023 
± 0.007 
± 0.004 
± 0.003 
2500 
± 0.008 



0.379 
79.29 
0.414 
0.0150 
7.700 
0.157 
0.211 
0.198 

0.870 



± 0.010 
± 0.16 
± 0.010 
± 0.0003 
± 0.019 
± 0.006 
± 0.003 
± 0.003 
2550 
± 0.008 



consistent within their statistical errors and a systematic 
uncertainty of the same order as the statistical one. Figure[8] 
shows the six different values of the volume-averaged sec- 
ondary star radius overplotted on a M-R relation for MS 
stars. Taken at face value, the results of the MCMC opti- 



misation indicate that the secondary is ~ 10 percent larger 
than theoretically predicted. As can be seen in Figure[8] this 
discrepancy drops to ~ 5 percent, if magnetic activity of 
the secondary is taken into account. With regard to the sec- 
ondary temperature, we note again that due to the black- 
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Figure 8. M ass-radius plot f or low -mass stars. Black point s 
are data from iLopez-Moraled l|2007t ) and iBeattv et all l|2007l) . 
where the masses of single stars were determined using mass- 
luminosity relations . The dotted line is the 5.0-Gyr isochrone 
from iBaraffe et al. I (QUI). The dashed li ne is a 5.0-Gyr mode l 
including effects of magnetic activity from Morales e t al] l|20ld ) . 
The results of the six chains for the volume-averaged radius of the 
secondary in SDSS1210 are plotted in red (online version only). 
Inset panel: zoom-in on the values of SDSS1210. The points are 
CIA: open circle; C1E: filled circle; C2A: open square; C2E: filled 
square; C3A: open triangle; C3E: filled triangle. 

body approximation, the value of T scc does not necessarily 
represent the true temperature of the star, it is effectively 
just a flux scaling factor. 

The gravitational redshift predicted by the light curve 
models (Table[4]), correcting for the redshift of the secondary 
star, the difference in transverse Doppler shifts and the po- 
tential at the secondary star owing to the white dwarf, 
are zwd = 15.9 ± 0.4 km s -1 from C1E and zwd = 

15.8 ± 0.4 km s -1 from C2E, where the errors are purely 
statistical and have been derived in the same manner as the 
other quantities reported in Tabled The systematic uncer- 
tainties in our photometric data might still be influencing 
the result, as the inclination angle and the stellar radii enter 
the calculation of zwd- Comparing zwd with the spectro- 
scopically determined value of zwd,s P cc = 7wd — 7scc = 

11.9 ± 1.7 km s" 1 we find that they are consistent within 
~ 2 <T. The systemic velocities 7wd and 7 sec are determined 
from spectroscopic observations obtained using a dual-arm 
spectrograph, with the white dwarf velocity measured in the 
blue arm and that of the secondary measured in the red arm 
(Sec.[2l[5]). The observations in both arms are independently 
wavelength-calibrated and the RMS of ~0.03A (Sec. [2]) cor- 
responds to an accuracy of the zero-point of ~ 1 — 2 kins -1 . 
The potential of an offset in the calibrations of the two arms 
enters the determination of zwd,s P cc as an additional sys- 
tematic uncertainty. 



8 PAST AND FUTURE EVOLUTION OF 
SDSS1210 

Considering its short orbital period, SDSS1210 must have 
form e d through com mon-envelope evolution JPaczvnski 
1 19761 ; IWebblnkll2008r i. As shown by ISchreiber fc Gansick" 



(2003), if the binary and stellar parameters are known, it 
is possible to reconstruct the past and predict the future 
evolution of PCEBs for a given angular momentum loss 
prescript ion. Here, we assume cl assical disrupted magnetic 
braking (|Verbunt fe Zwaanlll98ll ). In this context, given the 
low mass of the secondary, the only angular momentum loss 
mechanism for SDSS1210 is gravitational radiation. Based 
on the temperature and the mass of the white dwarf we in - 
terpolate the cooling tracks of lAlthaus fc Benvenutd (|l997T l 
and obtain a cooling age of t coo i = 3.5 Gyr. This corre- 
sponds to the time that passed since the binary left the 
common envelope. We calculate the period it had when 
it left the common envelope to be Pce = 4.24h. Fol- 
lowing the same method as in IZorotovic et~ai1 (|201Ch and 
based on their results we reconstructed the initial param- 
eters of the binary using a common-envelope efficiency of 
acE = 0-25 and the same fraction of recombination energy 
(see IZorotovic et alj 1201(3 . for more details). We found an 
initial mass of A/ prog = 1.33 Mq for the progenitor of the 
white dwarf, which filled its Roche lobe when its radius was 
^prog = 91.3 Rq. At that point, the orbital separation was 
a = 162.7 Rq, and the age of the system was t sys = 4.4 Gyr, 
since the time it was formed. Using the radius of the sec- 
ondarjjf] we calculate that the system will reach a semi- 
detached configuration and become a cataclysmic variable 
(CV) at an orbital period of P s( j ~ 2h in t st j = 1.5 Gyr. 

Given that the current P or b places SDSS1210 right 
at the upper edge of the CV orbital period garjf], and 
that the calculated P s< j, when SDSS1210 will start mass- 
transfer, is right at the lower edge of the period gap, we 
are tempted to speculate whether SDSS1210 is in fact a de- 
tached CV enterin g (or just having entered) the period gap. 
I Davis et al.l l|2008h have shown that a large number of de- 
tached WD+MS binaries with orbital periods between 2-3 
hours are in fact CVs that have switched off mass-transfer 
and are crossing the period gap. This could in principle ex- 
plain the apparently over-sized secondary in SDSS1210, as 
expected from the disr upted magnetic braking theory (e.g. 
iRappaport et all Il983l ). However, the temperature of the 
WD in SDSS1210 seems to be uncom fortably low for a WD 
that has recently stopped accreting (|Townslev Sz Gansickel 
120091 ). 



9 DISCUSSION AND CONCLUSIONS 

In this paper, we have identified SDSS1210 as an eclipsing 
PCEB containing a very cool, low-mass, DAZ white dwarf 
and a low-mass main-sequence companion. 

Using combined constraints from spectroscopic and 
photometric observations we have managed to measure the 
fundamental stellar parameters of the binary components. 
Systematic uncertainties in the absolute calibration of our 
photometric data, influence the determination of the stellar 
radii. The stellar masses, however, remain unaffected and 
were measured to a 1% accuracy. The (formal) statistical 
uncertainties in all binary parameters indicate the level of 

5 We assume a representative value of i? scc vol. aver. = 0-2 R© for 
the volume- averaged radius of the secondary 

6 The orbital period range where only a small number of CVs are 
found. 
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Table 5. Adopted stellar and binary parameters for SDSS1210. 



Parameter 



Value 



Porb [d] 

q 

o[Rq] 
Inclination [°] 
M WD [M Q ] 

%d [R©] 
log 9 

^eff.WD [K] 

if WD [kms -1 ] 

M scc [M ] 
-Rscc [R©] 

sec, vol. aver. [R©] 
Jf S ec[kms 1 ] 



0.124 489 764(1) 
0.379 ± 0.009 
0.870 ± 0.008 
(79.05 - 79.36) ± 0.15 
0.415 ± 0.010 
(0.0157 - 0.0161) ± 0.0003 
7.65 ± 0.02 
6000 ± 200 
95.3 ± 2.1 
0.158 ± 0.006 
(0.210 - 0.217) ± 0.003 
(0.197 - 0.203) ± 0.003 
251.7 ± 2.0 



precision that can be achieved in this system. All parameters 
are summarised in Table[5] 

With a mass of Mwd = 0.415 ± 0.010 M Q and a tem- 
perature of T e ff,WD ~ 6000 K, the DAZ white dwarf in 
SDSS1210 pushes the boundaries in a hitherto unexplored 
region of the WD parameter space. The M-R results from the 
four Chains CI and C2 are consistent with a He-core WD, 
assuming a hydrogen layer of M (H) /M WD = 3 x 10" 4 . How- 
ever, due to lack of observational constraints for the H-layer 
thickness and the uncertainty in the radii, we will defer iden- 
tifying the WD as a definite He-core and simply emphasise 
the strong candidacy. 

The secondary star, with a mass of M BCC = 0.158 ± 
0.006 M©, illustrates once more the excellent opportunity 
that PCEBs give us for testing and calibrating the M- 
R relations of low-mass stars. Taking the radius mea- 
surements at face value, the secondary star seems to be 
~ 10 percent larger than the theoretical values, although 
this drops to ~ 5 percent, if magnetic activity is taken 
into consideration. In this context, the magnetic activ- 
ity present in the secondary can lead to the formation 
of stellar (dark) spots on the surface. The effect of these 
spots is to block the outgoing heat flux, reducing T e fr 
and, as a result, the secondary expands to maintain ther - 
mal equilibr i um Ichabrier et al.ll2007l ; iMorales et~afll201Ch . 
iKraus et alj (|2011 ) found that low-mass stars in short pe- 
riod binaries appear to be overinflated (although their anal- 
ysis was restricted to M scc > 0.3 Mq), which seems to be 
the case for SDSS1210. We should note however, that the 
mass and radius of the secondary star in the eclipsing PCEB 
NNSer (with M BCC = 0.111 ± 0.004 M and comparable or- 
bital period to SDSS1210) is consistent with theoretical M-R 
predictions, e ven though it is heav ily irradiated by the hot 
WD primary (jParsons et al.ll2010al ). 

We have speculated whether SDSS1210 is in fact a de- 
tached CV entering the period gap, which could explain 
the large radius of the secondary. This hypothesis could 
be tested by measuring the rotational velocity of the white 
dwarf. This can be achieved through high-resolution spec- 
troscopy of_the_^e^_absorption lines in the WD photosphere 
see e.g. iTappert et afl|201ll ). 

In any case, it is highly desirable to improve the mea- 
surement of the stellar radii in SDSS1210 to the compara- 
ble precision to the masses presented here. This will require 



high-precision photometry in standard filters, such as e.g. 
delivered by ULTRACAM (|Dhillon et alj|2007h . 
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